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摘要 :【 目 的 】 本 研究 旨 在 建立 携带 与 未 携带 松 材 线虫 Bursaphelenchus xylophilus 的 松 墨 天 牛 
Monochamus alternatus 成 来 各 组 织 的 转录 组 数据 库 , 揭 示 松 墨 天 牛 对 松 材 线 虫 响 应 的 转录 组 整体 表 
达 特 征 。【 方 法 】 以 携带 和 未 携带 松 材 线虫 的 松 墨 天 牛 成 下 的 表皮 、 气 管 和 脂肪 体 为 材料 ,采用 
Illumina HiSeq™ 2000 测序 平台 开展 转录 组 测序 ,利用 Trinity 软件 对 RNA-Seq 数据 进行 从 头 组 装 ， 
利用 NCBI 数据 库 进行 基因 注释 。 通 过 DEG-Seq 分 析 携 带 和 未 携带 松 材 线虫 的 松 墨 天 牛 成 正中 的 
差异 表达 基因 ,对 上 调 表达 基因 进行 GO 和 KEGG 代谢 途径 富 集 分 析 , 并 通过 实时 荧光 定量 PCR 
技术 对 部 分 差异 表达 的 基因 包括 ENV, CDKI, HSP70-C, HSP70, HSP75, DUO, PABPNI 和 IGFP 
基因 的 表达 水 平 进行 验证 。【 结果】 测序 过 滤 后 共 获 得 55 059 条 单 基 因 往 (unigene) ,平均 长 度 为 
1 536 bpo 4% unigene 与 数据 库 中 的 序列 进行 BLASTX 比 对 ,成 功 注释 24 354 条 unigenes ,其 中 4 022 
条 unigene 在 GO 数据 库 中 获得 注释 ,6 098 条 unigene Æ KEGG 数据 库 中 获得 注释 。 经 过 G0 与 
KEGG 富 集 分 析 显 示 , 松 墨 天 牛 成 虫 携 带 松 材 线虫 后 ,差异 表达 基因 主要 为 与 其 气管 中 压力 调节 、 
组 织 与 DNA 修复 激素 响应 方面 相关 的 基因 。 实 时 荧光 定量 PCR 分 析 结 果 显 示 CDK1 ，HSP70， 
HSP75, DUO, PABPNI 和 IGFP 基因 的 表达 量 在 携带 松 材 线虫 的 松 墨 天 牛气 管 中 相 比 未 携带 松 材 
线 吕 的 松 墨 天 牛气 管 中 的 显著 上 调 。 【结论 】 本 研究 初步 阐明 松 墨 天 牛 携带 松 材 线虫 后 转录 组 的 
整体 表达 模式 ,为 进一步 研究 松 置 天 牛 的 基因 功能 及 对 松 材 线虫 的 胁迫 响应 过 程 黄 定 了 基础 。 
关键 词 : 松 墨 天 牛 ; 松 材 线虫 ; 转录 组 ; 基因 注释 ; GO; KEGG 
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Abstract: [ Aim] The objective of this study is to establish the transcriptome database of adult tissues of 
the pine sawyer beetle, Monochamus alternatus , carrying and without carrying the pine wood nematode, 
Bursaphelenchus xylophilus, and to reveal the overall expression characteristics of M. alternatus in 
response to pine wood nematode. [Methods] The transcriptomes of epidermis, trachea and fat body of 


M. alternatus adults carrying and without carrying B. xylophilus were sequenced using the Illumina 








基金 项 目 : 中 国 科学 院 前 沿 科学 重点 研究 项 目 (QYZDB-SSW-SMC014) ; 国家 自然 科学 基金 项 目 (31572272) 

作者 简介 : FBF, B, 1993 4E 2 HAE, 河北 沧州 人 , 硕士 研究 生 , 人 研究 方向 为 林业 有 害 生物 综合 防治 ,E-mail: 0713lmy@ gmail. com 
“通讯 作者 Corresponding authors, E-mail; whx25@ 126. com; zhaoll@ ioz. ac. cn 

收 稿 日 期 Received: 2019-09-02; 接受 日 期 Accepted: 2019-10-16 
































208 


昆虫 学 报 Acta. Entomologica Sinica 


63 卷 





HiSeq™ 2000 platform, the transcriptome data were de novo assembled using Trinity software, and gene 
annotation was conducted by using the NCBI database. The differentially expressed genes in M. 
alternatus adults carrying and without carrying B. xylophilus were screened by DEG-Seq, and the up- 
regulated genes were subjected to GO and KEGG metabolic pathway enrichment analysis. The expression 
levels of some differentially expressed genes including ENV, CDK1, HSP70-C, HSP70, HSP75, DUO, 
PABPNI and IGFP genes were verified by real-time quantitative PCR ( RT-qPCR). [Results] After 
sequencing and filtering, a total of 55 059 unigenes were obtained with an average length of 1 536 bp. A 
total of 24 354 unigenes were annotated using BLASTX searches against the databases, of which 4 022 
unigenes were assigned in GO database and 6 098 unigenes were assigned in the KEGG database. The 
enrichment analysis of GO and KEGG showed that the main differentially expressed genes in the trachea 
of M. alternatus adults carrying B. xylophilus are related with pressure regulation, tissue and DNA 
repair, and hormone response. The RT-qPCR analysis showed that the expression levels of CDKI, 
HSP70, HSP75, DUO, PABPNI, and ICFP genes were significantly higher in the trachea of M. 
alternatus carrying pine wood nematodes than in the trachea of M. alternatus without carrying pine wood 
nematodes. [ Conclusion] This study preliminarily clarified the overall expression pattern of the 
transcriptome of M. alternatus carrying B. xylophilus, providing a foundation for further study of the gene 
function of M. alternatus and the stress response of M. alternatus to B. xylophilus. 


Key words: Monochamus alternatus; Bursaphelenchus xylophilus ; transcriptome; gene annotation; GO; KEGG 





松 材 线虫 Bursaphelenchus xylophilus 属于 侧 尾 腺 
5j] ( Secernentea ) YẸ 7J H ( Aphelenchida ) i8 7] f] 
( Aphelenchoididae ) 4d um Bursaphelenchus 
( Nickle et al., 1981) ,是 一 种 林业 重大 外 来 入 侵害 
虫 ,造成 的 松 材 线虫 病 , 即 松树 的 “癌症 ”, 具 有 传染 
快 致死 快 ,防治 困难 ,一 经 感染 便 无 法 治愈 等 特点 
(Mamiya, 1972, 1983; THESE, 2005) 。 松 材 线虫 
原 产 于 北美 洲 ,于 20 世纪 初 传 人 亚洲 , 近 些 年 来 又 
扩散 到 了 欧洲 国家 , 除 北美 洲 原 发 地 以 外 ,对 许多 
家 的 森林 生态 系统 造成 了 严重 的 破坏 。 松 材 线 虫 于 
1982 年 首次 在 我 国 南 京 中 山陵 发 现 ,而 后 不 断 扩 
散 ,快速 蔓延 至 江苏 、 浙 江 、 安 徽 . 广 州 和 山东 等 53 
个 县 市 区 ,累计 死亡 松树 3 500 万 株 ,近年 来 疫情 逐 
步 加 剧 ,向 中 国 北部 划 延 (Yano et al., 1913; 孙 永 
Æ, 1982; Futai, 2008; Nakamura-Matori et al., 
2008; Shin, 2008; and Sahajdak, 
2010) 。 目 前 , 松 材 线虫 已 成 为 我 国 第 一 大 林业 害 
虫 ,也 是 我 国 最 为 严重 的 林业 外 来 有 害 生 物 ,给 病害 
REKER TIE WHR OKE EMKA K, 
2003). 。 随 着 国际 贸易 的 全 球 化 及 我 国 加 入 世贸 组 
织 , 松 材 线虫 由 其 他 疫 区 传人 我 国 非 疫 区 的 可 能 性 
不 断 加 大 ,风险 不 断 加 剧 ,而 且 由 于 国内 各 种 基建 项 
目的 施行 及 国内 经 济 贸易 活动 的 活跃 , 松 材 线虫 由 
我 国 国内 疫 区 传人 非 疫 区 的 可 能 性 也 在 逐渐 加 大 。 
因此 , 松 材 线虫 在 我 国 仍 具 有 十 分 广泛 的 定 殖 及 扩 





Karnkowski 











散 的 可 能 ( 谢 两 炎 等 , 2009; TRIER, 2015), 

在 中 国 和 日 本 , 松 材 线虫 的 感染 周期 与 松 墨 天 
^F. Monchamus alternatus 的 羽化 过 程 同步 ,并 且 没 有 
松 墨 天 牛 的 携带 , 松 材 线虫 便 无 法 被 传播 。 松 材 线 
虫 与 松 墨 天 牛 间 存 在 特殊 的 共生 生活 史 。 松 材 线虫 
生活 史 包 含 繁 殖 和 扩散 两 个 周期 ,繁殖 周期 主要 扩 
大 种 群 密 度 ,致死 松树 ;扩散 周期 则 主要 通过 媒介 昆 
虫 转移 至 新 寄主 。 当 松木 内 松 材 线虫 大 量 繁殖 ,其 
对 松树 的 危害 加 剧 ,造成 树 势 衰弱 甚至 枯死 ,从 而 吸 
引 松 墨 天 牛 前 来 产 卵 。 秋 末 冬 初 , 此 时 树 势 已 不 利 
于 繁殖 型 松 材 线虫 生长 ,繁殖 型 2 龄 幼虫 向 扩散 型 
3 龄 幼虫 (Ln ) 转变 。 此 时 , 松 墨 天 牛 孵 化 、 幼 虫 蚁 
皮 ; 翌 年 早春 , 树 体 内 松 墨 天 牛 进入 晴 期 , 松 材 线虫 
扩散 型 3 龄 幼虫 (Ln ) 向 松 墨 天 牛 晴 室 聚集 ,松树 蒸 
腾 作 用 降低 、 针 叶 凌 芒 黄 化 。 松 材 线虫 扩散 型 dt 
幼虫 (Lrm ) 回 扩散 型 4 龄 幼虫 (Lv ) 的 转变 与 松 墨 天 
牛 晴 的 发 育 期 相 一 致 ,翌年 春天 , 待 松 墨 天 牛 从 死 树 
中 羽化 ,可 耐 干旱 及 恶劣 环境 ,也 适合 通过 媒介 昆虫 
传播 的 松 材 线虫 扩散 型 4 龄 幼虫 (Lv ) 同步 形成 ,并 
通过 松 墨 天 牛气 管 转移 至 新 的 健康 松树 ,每 头 松 墨 
天 牛 平均 携带 量 约 18 000 条 ,最 多 可 达 289 000 条 
(Mamiya, 1983) 。 

松 材 线虫 与 松 墨 天 牛 的 生活 史 显 示 ,两 者 间 特 
殊 的 共生 关系 构成 了 松 材 线虫 病 形成 的 生物 学 基 
础 。 松 材 线虫 扩散 型 4 龄 幼虫 (Lv ) 通过 其 关键 媒 
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介 松 墨 天 牛 的 气管 被 携带 传播 ,而 松 墨 天 牛 与 松 材 
线虫 转型 扩散 发 育 进程 的 一 致 性 是 松 材 线虫 扩散 表 
型 增强 的 关键 。 松 墨 天 牛 在 羽化 过 程 中 释放 大 量 
Cis - Cis 脂肪 酸 乙 酯 ,启动 松 材 线虫 Insulin 和 DAZ 
DAF-12 信号 通路 ,改变 发 育 途 径 , 进 而 诱导 扩散 型 
3 龄 幼虫 (Ly ) 形 成 扩散 型 4 龄 幼虫 (Ly ) (Zhao et 
al., 2013) 。 由 Ly ZH DOSE RE UE US SEA SS D 
Bd e CR ^E XO E: Va] do Re t RR OS DSL U^ , 进 
而 促进 松 墨 天 牛 肾 化 。 一 旦 天 牛 进入 成 虫 阶段 ,天 
牛 分 泌 映 息 并 吸引 Liy ,从 而 促进 松 材 线虫 进入 松 墨 
天 牛气 管 并 被 转运 至 新 松树 上 (Zhao et al., 2016) 。 

总 之 ,深入 研究 松 墨 天 牛 与 松 材 线虫 “协同 互 
作 ” 关 系 , 对 解释 松 材 线虫 扩散 型 4 龄 幼虫 (Ly ) 被 
松 墨 天 牛 携带 机 制 具 有 重要 意义 ,对 控制 松 墨 天 牛 
携带 松 材 线虫 的 传播 具有 积极 的 理论 和 实际 意义 。 
因此 ,我 们 利用 高 通 量 umina HiSeq™ 2000 测序 技 
术 、. 通 过 组 装 .功能 注释 和 分 析 , 全 面 探讨 松 墨 天 牛 
携带 松 材 线虫 后 各 个 组 织 的 转录 组 信息 ,挖掘 松 墨 
天 牛 应 对 松 材 线虫 的 压力 响应 基因 ,从 分 子 水 平 上 
提供 松 墨 天 牛 携带 松 材 线虫 后 维持 稳 态 的 基因 信 
E ,也 为 深入 研究 松 墨 天 牛 能 够 携带 松 材 线虫 的 机 
制 英 定 基础 。 
1 材料 与 方法 
1.1 供 试 昆虫 

携带 松 材 线虫 的 松 黑 天 牛 采集 自 浙江 富阳 (30° 
07'N, 119°95 往 ) 死 亡 的 松树 上 ,对 照 组 所 用 的 未 携 
带 任何 线虫 的 松 墨 天 牛 收集 自 实 验 室 种 群 。 实 验 室 
人 工 饲养 条 件 : 幼 虫 放 在 装 有 木 导 饲料 的 10 cm 的 
指 型 管 中 , 每 7 d 换 一 次 饲料 ;成 虫 置 于 塑料 盒 中 ， 
每 7 d 换 一 次 新 鲜 松 枝 , 以 上 培养 过 程 均 置 于 人 工 
气候 箱 中 ,饲养 温度 为 20 ~ 25Y ,相对 湿度 60% ~ 
7096 , 光 周 期 12L: 12D。 所 有 天 牛 在 实验 前 均 解剖 
出 气管 并 检测 携带 线虫 情况 及 数量 后 ,再 收集 组 织 
用 于 文库 构建 。 所 选 携带 松 材 线虫 的 松 墨 天 牛 样 
本 ,体内 携带 松 材 线虫 的 数量 为 1 000 ~3 000 头 。 
1.2 文库 构建 

收集 未 携带 松 材 线虫 的 松 墨 天 牛 成 虫 表皮 
(EP-CK) .气管 (Tr-CK) .脂肪 体 (Fb-CK) 以 及 携带 
松 材 线虫 的 松 墨 天 牛 表皮 (EP-Bx) 脂肪 体 (Fb-Bx) 
和 气管 (Tr-Bx), 用 PBS DE TK, 随后 使 用 Trizol 
(Thermo Fisher Scientific ) 从 组 织 中 分 离 总 RNA ,并 
在 使 用 前 储存 在 - 80% 冰箱 。 共 构建 6 个 基因 文 
库 , 每 个 文库 构建 使 用 5 头 天 牛 的 组 织 合 并 为 一 个 























































































































样品 。 剩 余 的 天 牛 尸 体温 入 PBS 中 统计 线虫 的 数 
量 。 用 ND-1000 分 光 光 度 计 (NanoDrop 
Technologies, Wilmington, DE, 美国) 测定 RNA 样 
品 浓度 和 纯度 。 用 Dynabeads mRNA 纯化 试剂 盒 
( Ambion; Thermo Fisher Scientific, 美国 ) 纯化 
mRNA, RNA-Seq 文库 制备 试剂 盒 ( Compatible with 
Illumina, San Diego, CA, 美国, Ilumina Whole- 
Genome Gene Expression BeadChip ) 用 于 转录 组 文库 
FJ ££, (E FH TE VETE R2 ( dT) R Ik ( Thermo Fisher 
Scientific, 美国 ) ARE Poly A 的 mRNA ,然后 对 其 
进行 片段 化 , 反 转 录 合 成 cDNA。 使 用 MiniElute XE 
胶 提 取 纯 化 试剂 盒 320 - 420 bp 产物 (Qiagen， 
Hilden, 德国 ) ,测序 是 在 中 国 科学 院 基 因 组 研究 所 
的 Ilumina HiSeq 2000 测序 平台 上 完成 。 
1.3 转录 组 拼接 及 注释 

在 组 装 之 前 ,去 除 低 质量 的 原始 reads, 过滤 去 
除 微生物 (细菌 /病毒 /真菌 ) 和 线虫 污染 ,获得 的 
clean reads 用 于 后 续 拼 接 。 短 reads 在 Trinity ( v. 
2013) 中 进行 De novo 组 装 ( Grabherr et al., 2011) 。 
组 装 的 序列 进一步 通过 CD-HIT 软件 (Li and 
Godzik, 2006) 进行 聚 类 ,以 精简 元 余 序 列 并 人 工 校 
对 。 为 了 使 转录 本 尽 可 能 全 面 并 避免 增加 组 装 的 复 
杂 性 ,将 6 个 库 合 并 组 装 , 使 用 Burrows-Wheeler 
Aligner( BWA) 软件 (Li and Durbin, 2009) 将 测序 读 
数 重新 整合 到 unigenes 上 ,评估 组 装 序列 的 准确 性 。 

使 用 BlastX 软件 ,将 拼接 转录 本 与 NCBI( http: 
// www. ncbi. nlm. nih. gov/) JE 26 4x 4& A c JE Fe 
(NR) 和 Swiss-Prot 数据 做 比 对 ,设置 E-value «10 7 
为 阔 值 ,以 第 一 匹配 序列 的 信息 作为 注释 结果 。 使 
用 Trinity 软件 中 的 RSEM 评估 基因 丰 度 水 平 , 同 时 
使 用 具有 P «0.001 的 DEG-Seq 软件 鉴定 差异 表达 
的 基因 (DEC ) ,并 使 用 以 下 规则 标记 。 将 FDR 
(false discovery rate ) 作为 筛选 的 关键 指标 ,将 1log， 
Ratiol Z2 H. Q-value x0. 05 的 基因 作为 差异 表达 基 
,计算 出 上 调 与 下 调 的 基因 数量 。 用 KEGG( http: 
// www. kegg. jp) 和 GO ( http: // wego. genomics. org. 
cn/ ) 数据库 ( Ogata et al., 1999) 进行 功能 富 集 分 析 。 
使 用 R 软件 Pheatmap 包 ( http: // www. r-project. 
org/; R Foundation for Statistical Computing, Wien, 
澳大利亚 ) 绘 制 热 图 。 本 研究 测 得 的 转录 组 数据 已 
上 传 美国 国家 生物 技术 信息 中 心 (NCBI) SRA 数据 
库 ,BioProject PRJNA374773 。 
1.4 实时 荧光 定量 PCR 验证 差异 表达 基因 的 表达 

采用 实时 荧光 定量 PCR 的 方法 检测 1.3 45 $8 












































210 昆虫 学 报 Acta. Entomologica. Sinica 


63 卷 





选 的 在 携带 与 未 携带 松 材 线虫 的 松 墨 天 牛 间 差 异 表 
达 的 基因 在 松 墨 天 牛 不 同 组 织 中 的 相对 表达 量 。 这 
些 基因 包括 气管 中 表达 上 调 的 核酸 内 切 酶 V 
(endonuclease V) 基因 (ENTY) 细胞 周期 蛋白 依赖 性 
激酶 1( cyclin-dependent kinase 1) 基因 ( CDK1)、 热 
激 和 蛋白 70-C( heat shock protein 70-C ) 基因 (HSP70- 
C) 、 热 激 和 蛋白 70 (heat shock protein 70) 基 
(HSP70) 双重 氧化 酶 (dual oxidase ) 基因 (DUO) 、 
多 腺 苷 酸 结 合 蛋 白 (PABPNI protein) 基 
(PABPN1) 胰岛 素 样 生 长 因子 结合 蛋白 复合 物 酸 
不 稳定 亚 基 (insulin-like growth factor-binding protein 
complex acid labile subunit) 基因 (JGFP)、 热 激 和 蛋白 
75 (heat shock protein 75) 基因 (HSP75 ) 等 。 使 用 
Primer Premier 6 设计 内 参 基 因 Actin 引物 序列 和 目 
的 基因 引物 序列 ,以 1.2 节 合 成 的 cDNA 为 模板 , 引 
物 信息 详 见 表 1。 所 需 引 物 由 北京 六 合 华 大 基因 科 
技 有 限 公 司 合成 。PCR 反应 体系 : 2 x Talent qPCR 
Premix 10 kuL， 上 下 游 引物 (10 hmolL) 各 0.8 uL, 
cDNA 2 pL, 超 纯 水 补 至 20 pL。 反 应 在 ABI 
7300Plus 荧光 定量 PCR 仪 上 进行 。 反 应 程序 : 95%C 
预 变 性 30 s; 95C 变性 5 s, 55'C 退火 30 s, 72%C 延 
伸 30 s, 40 次 循环 ;接着 进行 溶解 曲线 分 析 。 表 达 
结果 依据 Ct 值 , 采用 2-^^" 法 计算 相对 表达 量 。 
使 用 SPSS19.0 软件 对 所 得 数据 进行 统计 分 析 。 









































2 结果 


2.1 RNA-Seq 数据 质 控 与 评估 

未 携带 松 材 线虫 的 松 黑 天 牛气 管 、 表 皮 及 脂肪 
体 的 Ilumina 测序 共 得 到 60 734 045 条 raw reads, 
去 除 不 可 用 数据 509 836 条 ,最 后 得 到 60 224 209 








条 clean reads。 携 带 松 材 线虫 的 松 墨 天 牛气 管 、 表 
皮 及 脂肪 体 的 Tlumina 测序 共 得 到 66 296 151 条 
raw reads, 去 除 不 可 用 数据 715 587 条 ,最 后 得 到 
65 580 564 条 clean reads( 表 2)。 上 述 结果 说 明 可 
用 reads 在 9996 以 上 ,本 研究 测 得 的 转录 组 数据 质 
量 优 ,可 用 于 进一步 分 析 。 

对 转录 本 进行 聚 类 和 组 装 分 析 共 得 到 55 059 
条 unigenes ,最 大 长 度 2 845 bp ,最 小 长 度 为 300 bp, 
平均 长 度 为 1 536 pp ,中 位 数 长 度 为 906 bp, N50 长 
度 为 2 516 bp。 超 过 1 kb 的 unigenes 有 12 394 条 ， 
占 22. 5% ,超过 2 kb 的 unigenes 有 13 219 条 , 占 
24.096 。 这 些 数 据说 明 组 装 效果 理想 。 











表 1 引物 信息 
Table 1 Primer information 
引物 引物 序列 (5' -3^) 

Primers Primer sequences 
ENV-F CAGAGCAGAACAACATATTAGG 
ENV-R GAGTTATTCCACAGACAGACT 
CDKI-F TGGCAATATCTCTCACGAA 
CDKI-R CGCCTTCATCTTCTTGTTC 
HSP 70-C-F CGTTCAAGGTGGAGCTTATCT 
HSP70-C-R GGTTATCAGCGGCAGTAG 
HSP 70-F AGTGTTGACCTGAAGAGAT 
HSP70-R GCTGTTCTTGAAGTGCATT 
DUO-F ATGACACTGCCTTACAAGAT 
DUO-R ATCGTTATGCTGCTCAACT 
PABPNI-F TCTTCGTTCTCTTCTCTGCTT 
PABPNI-R GAGGATGTGGCTGCTATTCA 
IGFP-F ACGAGGCAATAATGTGGAAGG 
IGFP-R GAACTGACGATAAGCAATAGC 
HSP75-F TCCTGCTCAGAGTACAATG 
HSP75-R CGTTCCTCGGCAAGATTA 
Actin-F TGGGTATGCAATCTTGCGGT 
Actin-R GGCCCTGATTTCCTTTTGCA 


表 2， 松 墨 天 牛 成 虫 各 组 织 转录 组 数据 评估 


Table 2 Analysis of the transcriptome data of adult tissues of Monochamus alternatus 























样品 原始 reads 数 低 质量 数据 数 微生物 (病毒 /细菌 ) 污染 数据 数 有 效 reads 数 有 效 reads 所 占 比例 (% ) 
Sample Number of Number of adaptors/ Number of reads of microbes Number of Proportion of 
raw reads low-quality reads ( viruses/ bacteria ) valid reads valid reads 

Ep-CK 22 359 506 190 048 8 965 22 160 493 99.11 

Tr-CK 19 190 469 171 267 7 938 19 011 264 99.07 

Fb-CK 19 184 070 123 780 7 838 19 052 452 99.31 

Ep-Bx 24 541 598 150 433 26 348 20 288 739 99. 11 

Tr-Bx 21 283 674 146 093 225 260 20 912 321 98.26 

Fb-Bx 20 470 879 155 792 11 661 24 379 504 99.34 








Ep-Bx; 携带 松 材 线虫 的 松 墨 天 牛 成 虫 表皮 Epidermis of M. alternatus adults carrying B. xylophilus; EP-CK: 未 携带 松 材 线虫 的 松 墨 天 牛 成 虫 表 























Hz Epidermis of M. alternatus adults without carrying B. xylophilus; Tr-Bx; 携带 松 材 线 虫 的 松 墨 天 和 牛 成 虫 气 管 Trachea of M. alternatus adults 














carrying B. xylophilus; Tr-CK: 未 携带 松 材 线虫 的 松 墨 天 牛 成 虫 气 管 Trachea of M. alternatus adults without carrying B. xylophilus; Fb-Bx: 携带 松 
材 线虫 的 松 墨 天 牛 成 虫 脂 肪 体 Fat body of M. alternatus adults carrying B. xylophilus; Fb-CK: 未 携带 松 材 线虫 的 松 墨 天 牛 成 虫 脂 肪 体 Fat body of 
M. alternatus adults without carrying B. xylophilus. 图 5 -6 和 表 3 [i] The same for Figs. 5 -6 and Table 3. 
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2.2 携带 和 未 携带 松 材 线虫 的 松 墨 天 牛 成 虫 组 织 
转录 组 中 的 差异 表达 基因 

以 未 携带 松 材 线虫 的 松 墨 天 牛 成 虫 各 组 织 的 转 
录 组 数据 为 对 照 ,携带 松 材 线虫 天 牛 成 虫 的 气管 中 
有 差异 表达 基因 7 089 个 ,其 中 上 调 基 因 4 526 个 ， 
下 调 基 因 2 563 个 。 表 皮 中 有 差异 表达 基因 3 923 
个 ,其 中 上 调 基 因 2 575 个 ,下 调 基 因 1 348 个 。 脂 
肪 体 中 有 差异 表达 基因 5 485 个 ,其 中 上 调 基因 
1 748 个 ,下 调 基 因 3 737 个 。 气 管 与 表皮 中 共有 的 
差异 表达 基因 1 871 个 ;其 中 二 者 全 部 上 调 的 基因 
849 个 ,全 部 下 调 的 基因 425 个 ;气管 中 上 调 而 表皮 
中 下 调 的 基因 265 个 ,气管 中 下 调 而 表皮 中 上 调 的 
基因 332 个 。 表 皮 和 脂肪 体 中 共有 的 差异 表达 基因 
1 587 个 ;其 中 二 者 全 部 上 调 的 基因 480 个 ,全 部 下 
调 的 基因 676 个 ;表皮 中 上 调 而 脂肪 体 中 下 调 的 基 
214 个 ,表皮 中 下 调 而 脂肪 体 中 上 调 的 基因 217 
个 。 气 管 与 脂肪 体 中 共有 的 差异 表达 基因 1 683 
个 ;其 中 二 者 全 部 上 调 的 基因 374 个 ,全 部 下 调 的 基 
545 个 ;气管 中 上 调 而 脂肪 体 中 下 调 的 基因 481 
个 ,气管 中 下 调 而 脂肪 体 中 上 调 的 基因 283 个 。 气 
E 表皮 与 脂肪 体 中 共有 的 差异 表达 基因 1 446 个 ; 
三 者 全 部 上 调 的 基因 为 488 个 ,全 部 下 调 的 基因 为 
367 个 ;气管 与 表皮 中 上 调 而 脂肪 体 中 下 调 的 基因 
为 103 个 ,气管 与 脂肪 体 中 上 调 表达 而 表皮 中 下 调 
表达 的 基因 为 78 个 ;表皮 与 脂肪 体 中 上 调 而 气管 中 
下 调 的 基因 为 89 个 ;气管 中 上 调 而 表皮 和 脂肪 体 中 
下 调 的 基因 为 163 个 ,表皮 中 上 调 而 气管 与 脂肪 体 
中 下 调 的 基因 为 58 个 ;脂肪 体 中 上 调 而 气管 与 表皮 
中 下 调 的 基因 为 100 个 (图 1)。 以 上 结果 说 明 , 松 
墨 天 牛气 管 相对 于 表皮 及 脂肪 体 是 主要 响应 松 材 线 
虫 的 组 织 。 

2.3 差异 表达 基因 的 功能 注释 

利用 BLASTX 程序 将 55 059 条 unigenes 序列 与 
NR(14 234 条 unigenes) , GO(4 022 条 unigenes ) 和 
KEGG(6 098 条 unigenes ) 数据 库 进行 比 对 ,成 功 注 
释 到 24 354 条 unigenes , E 44.2396 , 

KEGG 数据 库 注释 到 6 098 条 基因 , 占 11.196 。 
属于 气管 与 表皮 中 的 共有 2 696 条 基因 ;经 过 比 对 ， 
相 较 于 未 携带 松 材 线虫 的 松 墨 天 牛 ,携带 松 材 线虫 
的 松 墨 天 牛气 管 中 上 调 表达 基因 与 下 调 表达 基因 最 
多 的 均 为 复制 和 修复 (replication and repair) 相关 
的 ,分 别 为 114 和 152 条 ;表皮 组 织 中 上 调 表达 基因 
最 多 的 为 折 厂 、 分 类 和 降解 (folding，sorting and 
degradation) 相关 的 ,为 54 条 ,下 调 表达 基因 最 多 的 
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图 1 携带 和 未 携带 松 材 线虫 的 松 墨 天 牛 成 虫 
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wW^100 YAY 163 
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转录 组 差异 表达 基因 韦 恩 图 


Fig. 1 Venn diagram of differentially expressed genes in 





Monochamus alternatus adults carrying and without 
carrying Bursaphelenchus xylophilus 
Tr; 气管 Trachea; Ep; 表皮 Epidermis; Fb; 脂肪 体 Fat body. 向 上 和 
向 下 的 箭头 分 别 表 示 上 调和 下 调 基因 。Up and down arrows indicate 


up-regulated and down-regulated genes, respectively. 














为 信号 分 子 与 相互 作用 (signaling molecules and 
interaction ) 相关 的 ,为 44 条 (图 2)。 

GO 数据 库 注 释 到 4 022 条 基因 , 占 7.3% ; 相 较 
于 未 携带 松 材 线虫 的 松 墨 天 牛 ,携带 松 材 线虫 的 松 
墨 天 牛 表皮 与 气管 中 上 调 的 基因 (图 3) 中 属于 细胞 
组 分 (cellular component ) 分 类 的 主要 有 细胞 (cell) 、 
细胞 外 区 域 组 分 (cell part ) 和 细胞 内 区 域 组 分 
(intracellular part) 等 ; 属于 分 子 功能 (molecular 
function ) 分 类 的 主要 为 催化 活性 (catalytic activity) 、 
水 解 酶 活性 (hydrolase activity ) 和 转移 酶 活性 
( transferase activity ) 等 ;属于 生物 学 过 程 (biological 
process) 分 类 的 有 代谢 过 程 (metabolic process ) 、 含 
氮 化 合 物 代 谢 过 程 (nitrogen compound metabolic 
process) 和 有 机 物 代 谢 过 程 (organic substance 
metabolic process) 等 。 下 调 的 基因 (图 4) 中 属于 细 
胞 组 分 分 类 的 主要 有 细胞 (cell) 、 细 胞 外 区 域 组 分 
(cell part) 和 细胞 内 液 (intracellular fluid) 等 ;属于 分 
子 功 能 分 类 的 主要 为 催化 活性 (catalytic activity ) 、 
结合 ( binding ) 4I Zi 3$ 3 DS] fj £s 4r. ( heterocyclic 
compound binding) 等 ;属于 生物 学 过 程 分 类 的 有 代 
谢 过 程 (metabolic process ) 、 含 所 化 合 物 代谢 过 程 
(nitrogen compound metabolic process) 和 有 机 物 代 谢 


过 程 (organic substance metabolic process) 等 。 
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图 3 在 携带 松 材 线虫 松 墨 天 牛 成 虫 表皮 与 气管 转录 组 中 上 调 基因 的 GO 注释 分 布 


Fig. 3 GO annotation distribution of up-regulated genes in epidermis and trachea transcriptomes of Monochamus alternatus 








adults carrying Bursaphelenchus xylophilus 
Tr-u; 气管 中 上 调 表达 的 基因 Up-regulated genes in trachea; Ep-u: 表皮 中 上 调 表达 的 基因 Up-regulated genes in epidermis. 
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图 4 在 携带 松 材 线虫 松 墨 天 牛 成 虫 表皮 与 气管 转录 组 中 下 调 基因 的 GO 注释 分 布 
Fig. 4 GO annotation distribution of down-regulated genes in epidermis and trachea transcriptomes of Monochamus alternatus 
adults carrying Bursaphelenchus xylophilus 
Tr-d: 气管 中 下 调 表 达 的 基因 Down-regulated genes in trachea; Ep-d: 表皮 中 下 调 表达 的 基因 Down-regulated genes in epidermis. 








其 中 ,GO 功能 富 集 分 析 结 果 表 明 , 气 管 样品 中 
上 调 的 基因 主要 富 集 在 分 子 功能 中 的 催化 活性 、 结 
合生 物 学 过 程 中 的 代谢 过 程 和 细胞 过 程 (图 3) ;下 
调 的 基因 则 主要 富 集 在 细胞 组 分 中 的 细胞 、 细 胞 外 
区 域 组 分 细胞 内 区 域 组 分 以 及 分 子 功能 分 类 中 的 


结合 杂 环 基因 的 结合 、 环 类 有 机 物 的 结合 和 离子 作 
用 (图 4)。 表 皮 样 品 中 上 调 的 基因 同样 主要 定 集 在 
分 子 功 能 中 的 催化 活性 、 结 合 以 及 生物 学 过 程 中 的 
代谢 过 程 有 机 物 代谢 过 程 及 初级 代谢 过 程 (图 4); 
下 调 的 基因 则 富 集 在 分 子 功能 中 的 催化 活性 、 结 合 
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以 及 生物 过 程 中 的 代谢 过 程 、 含 氮 化 合 物 代谢 过 程 
和 有 机 物 代谢 过 程 (图 4)。 
2.4 松 墨 天 牛 胁迫 响应 相关 基因 的 鉴定 
iit GO, KEGG 和 NR 注释 分 析 显 示 , 松 墨 天 
牛 携带 松 材 线虫 后 ,复制 和 修复 相关 功能 通路 基因 
达 显著 上 调 , 热 图 分 析 进 一 步 表 明 其 在 携带 松 材 
线虫 的 松 墨 天 牛 的 表皮 和 和 气管 中 显著 上 调 ( 图 5)。 
其 中 ,气管 中 上 调 表 达 的 共 济 失调 蛋白 (ataxin-3- 
like) 基因 8:38 E A (zinc finger protein) 基因 、 小 核 
糖 核 蛋白 颗粒 蛋白 SmE (small ribonucleoprotein 
particle protein SmE ) 基因 , IGFP 基因 、ENV 基因 、 
CDK1 基因 、HSP70-C 基因 、HSP70 基因 、DUO 基因 、 
中 链 酰 基辅 酶 A 脱氧 酶 (medium-chain specific acyl- 
CoA dehydrogenase, mitochondri ) 基因 、 锚 蛋白 重复 和 
FYVE W %Æ AĤ (ankyrin repeat and FYVE domain- 
containing protein 1) 基因 、CREB 结合 蛋白 (CREB- 
binding protein) 基因 ,参与 者 泡 运输 的 内 质 网 (ER) 
蛋白 (vesicle-associated membrane protein-associated 
protein B) 基因 , HSP75 基因 ,可 能 的 RNA 解 旋 酶 
(putative RNA helicase mov-10-B. 1) Æ KRF MA 
HRE A REMIJA] ( Cystatin-C) 基 因 等 ,与 胁迫 响应 相 
关 ( 表 3)。 经 过 定量 PCR 验证 ,其 中 CDKI 基因 、 
HSP70 基因 、HSP75 基因 、DUO 基因 、PABPN1 基因 和 
IGFP 基 因 的 表达 量 在 携带 松 材 线虫 的 松 墨 天 牛 成 


























虫 气管 中 相 比 未 携带 松 材 线虫 的 松 墨 天 牛气 管 中 的 
HEIE 6)。 
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图 5 胁迫 响应 相关 基因 在 携带 与 未 携带 松 材 线虫 的 
松 墨 天 牛 成 虫 组 织 中 的 表达 热 图 


Fig. 5 Heat map showing expression of genes related 











to stress response in the adult tissues of Monochamus 
alternatus carrying and without carrying 


Bursaphelenchus xylophilus 


表 3 携带 和 未 携带 松 材 线虫 的 松 墨 天 牛 成 虫 各 组 织 中 显著 上 调 的 压力 调节 相关 基因 


Table 3 Significantly up-regulated genes related to pressure regulation in adult tissues of Monochamus alternatus 


carrying and without carrying Bursaphelenchus xylophilus 













































































基因 ID NR 数据 库 注 释 信 息 基因 表达 量 Gene expression level 
Cene ID Nr database annotation information Ep-Bx Tr-Bx Fb-Bx Ep-CK . Tr-CK Fb-CK 

comp66412 c0 seql3 ”双重 氧化 酶 DUOX 0.43 1.02 0.42 0.59 0.68 0 
胰岛 素 样 生 长 因子 结合 蛋白 复合 物 酸 不 稳定 亚 基 

comp71324 c0. seql Insulin-like growth factor-binding protein complex acid 1.95 2.47 0.68 1.37 0.22 0.59 
labile subunit 
组 蛋白 赖 氨 酸 N- 甲 基 转 和 

comp1416939 cO. seql am A ERE 0 0.79 0.16 0 0 0 
Histone-lysine N-methyltransferase 

comp56144 c0 seq4 核酸 内 切 酶 V Endonuclease V 0.11 2:2 0.09 0 0 0 
线粒体 中 链 酰基 辅酶 A 脱 氢 酶 

comp459307. c0. seql Medium-chain specific acyl-CoA dehydrogenase, 0.21 1.82 0 0 0 0 
mitochondrial 

嘎 吟 无 喀 啶 DNA 裂 

comp38311. cO, seql nius . T 裂解 本 0. 19 1.26 0 0 0 0 
DNA-( apurinic or apyrimidinic site) lyase 

comp473734 c0. seql DNA 连接 酶 1 DNA ligase 1 0.22 1.81 0 0 0 0 
A Gate S80 AURI CER) EH B 

comp61820_c0_seql SEMAR ER RAB | 0.56 8.54 0 0 0 0 
Vesicle-associated membrane protein-associated protein B 
Z gE; ^M EET 1 

comp557178. c0. seql MERRE 1 依赖 性 激酶 4.05 4.57 1.97 3.8 0.93 1.34 
Cyclin-dependent kinase 1 
3 ^ 75 kDa ,线粒体 样 

comp616611_c0_seql SH a RAA 0 2.21 0 0 0 0 


Heat shock protein 75 kDa，mitochondrial-like 












































图 6 胁迫 响应 相关 基 























Tr-CK Tr-Bx 
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续 表 3 Table 3 continued 
基因 ID NR 数据 库 注 释 信 息 基因 表达 量 Cene expression level 
Gene ID Nr database annotation information Ep-Bx Tr-Bx Fb-Bx Ep-CK Tr-CK Fb-CK 
一 般 转 录 和 DNA 修复 因子 IIH 解 旋 酶 亚 基 XPD 
comp2097377_c0_seql General transcription and DNA repair factor HH helicase 1.48 3.43 0.12 2.89 0.76 0.38 
subunit XPD 
comp677104 c0 seql 蛋白 质 IMPACT 同 源 物 Protein IMPACT homolog 0. 34 2.64 0 0 0 0 
comp57202. c0. seql 可 能 的 RNA 解 旋 酶 Putative RNA helicase mov-10-B. 1 0. 12 4.08 0.16 0 0 0 
comp58077. c0. seql 半 胱 氨 酸 蛋白酶 抑制 剂 Cystatin-C 0.21 6.49 0.25 0 0 0 
comp395508 c0 seql IRANIN 1 Flap endonuclease 1 0.06 1.32 0 0 0 0 
comp63149 cO seql 热 激 和 蛋白 40 Heat shock protein 40 1.1 15.54 0.28 0 0 0 
comp45847. c0. seql 多 腺 苷 酸 结合 蛋白 PABPNI protein 0.56 8.54 0 0 0 0 
小 核糖 颗粒 蛋白 SmE 
comp134747. c0. seql 核糖 核 入 | WEI l pi . 0.59 6.72 0.16 0 0 0 
Small ribonucleoprotein particle protein SmE 
A 6 B 10 
E 58 了 
i E 4 i E 
2i Sis 
Ro Ro 
> i. 
a= "34 
gS 2 E. 
PER: PER: 
名 8 5 
0 0 
Tr-CK Tr-Bx Tr-CK Tr-Bx 
C 25 Ak D 6 
A20 EB 
iN E iN E 4 
EE EH 
Ro Ro 
TE T 
5-210 pe um 
m.s m.s 2 
PER: E-ER: 
B 5 [zi 
0 0 
Tr-CK Tr-Bx Tr-CK Tr-Bx 
E 20 F 5 
ek i 
£15 Eos 
E E 
3 8 1 E3 
To 10 Sog 
E * E 42 
x . X 
[a [a 1 
0 0 
Tr-CK Tr-Bx Tr-CK Tr-Bx 
G 25 H g 
* 
z20 s EB 6 
d m ii g 
Ro Ros 
$. pcd 
5-210 a 5 
uu mu 
Ex E, 
B 5 名 
0 0 


Tr-CK Tr-Bx 





因 在 携带 与 未 携带 松 材 线虫 的 松 墨 天 牛 成 虫 气管 中 表达 的 实时 荧光 定量 PCR. 验证 


Fig. 6 Expression verification of stress response-related genes in trachea of Monochamus alternatus adults carrying and 


A: ENV; B: CDK1; C: HSP70-C; D; HSP70; E: DUO; F; PABPNl ; G: IGFP; H; HSP75. 基因 表达 量 以 Tr-CK( 对 照 ) 中 的 为 基准 , 柱 上 


without carrying Bursaphelenchus xylophilus by real-time quantitative PCR 

















和 双星 号 分 别 表 示 两 组 间 差 异 显著 (P<0.05) 和 极 显 著 (P <0.01)( 独 立 样本 :检验 )。The gene expression level in Tr-CK ( control) was 
normalized to 1. The asterisk and double asterisk above bars indicate significant difference ( P <0.05) and extremely significant difference (已 <0.01) ， 


respectively, between 





two groups (independent samples t-test ) . 
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3 讨论 


目前 ,在 全 球 范围 内 植物 病原 体 已 经 带 来 了 许 
多 毁灭 性 的 病害 ,而 昆虫 媒介 负责 传播 一 些 世 界 上 
最 致命 的 病原 体 , 且 某 些 病原 需要 特定 的 昆虫 才能 
传播 (Suzuki et al., 2006) 。 媒 介 昆 虫 对 植物 病原 体 
的 携带 主要 集中 在 果 晶 或 半 却 目 昆虫 的 肠 道 和 唾液 
腺 中 (Basset et al., 2000) 。 昆 虫 对 线虫 的 携带 研究 
主要 在 昆虫 对 寄生 性 线虫 和 昆虫 病原 线虫 之 间 展 开 
( Castillo et al., 2011) 。 然 而 ,媒介 昆虫 对 植物 病原 
线虫 的 携带 ,尤其 是 在 气管 系统 中 则 人 鲜 有 报道 。 媒 
介 昆 虫 松 墨 天 牛 能 安全 地 携带 成 千 上 万 头 植物 病原 
线虫 松 材 线虫 而 不 影响 其 寿命 ,所 以 研究 松 墨 
天 牛 与 松 材 线虫 间 的 共生 互 作 关系 十 分 必要 。 本 人 研 
究 利 用 Humina 测序 技术 对 携带 与 未 携带 松 材 线虫 
的 松 墨 天 牛 进 行 转录 组 测序 和 分 析 ,在 无 参考 基因 
组 转录 组 测序 情况 下 , 共 获 得 55 059 条 unigenes。 
其 中 4 022 条 unigene 在 GO 数据 库 中 获得 注释 ， 
6 098 条 unigene 在 KEGG 数据 库 中 获得 注释 。 通 过 
分 析 后 发 现 松 墨 天 牛 携带 松 材 线虫 后 ,差异 基因 主 
要 为 松 墨 天 牛气 管 中 压 力 调节 组织 与 DNA 修复 相 
关 通 路 ,为 进一步 发 据 重 要 的 功能 基因 奠定 基础 。 
通过 对 比 携带 及 未 携带 松 材 线虫 的 松 墨 天 牛 的 
气管 .表皮 及 脂肪 体 的 转录 组 数据 ,发 现 除 脂 肪 体 
外 ,携带 松 材 线虫 的 松 墨 天 牛气 管 与 表皮 组 织 中 的 
差异 表达 基因 均 以 上 调 为 主 ,其 中 尤 以 气管 中 的 基 
因 表 达 量 变化 最 大 ,高 于 其 余 两 组 织 中 的 差异 表达 
基因 总 数 及 上 调 基 因数 。 对 差异 表达 基因 的 GO 和 
KEGG 分 析 结 果 显 示 ,携带 松 材 线虫 的 松 墨 天 牛气 
管 及 表皮 中 富 集 到 的 功能 通路 相 一 致 , 且 表 达 模 式 
也 相 一 致 。 说 明 尽 管 松 材 线虫 在 天 牛气 管 中 被 携 
带 , 昆 虫 气 管 系统 是 由 外 胚层 细胞 内 隐形 成 ,其 内 壁 
与 表皮 相连 ,与 表皮 构造 相同 ,昆虫 表皮 及 气管 系统 
都 具有 一 层 上 皮 细 胞 且 相 连 ,作为 昆虫 的 第 一 道 防 
线 抵挡 外 来 物 侵 染 。 前 期 研究 显示 , 松 墨 天 牛 对 松 
材 线虫 的 响应 过 程 中 ,总 体 上 松 黑 天 牛 的 免疫 处 于 
压制 状态 ,表皮 和 气管 中 一 部 分 免疫 基因 上 调 表 达 ， 
大 部 分 免疫 基因 表达 水 平 不 变 , 松 墨 天 牛 通过 活性 氧 
调节 维持 二 者 间 免 疫 稳 态 (Zhou et al., 2017, 2018), 
在 本 研究 中 ,我 们 通过 进一步 转录 组 分 析 ,筛选 到 大 
量 胁 迫 啊 应 基因 在 携带 松 材 线虫 的 松 墨 天 牛 中 上 
调 , 说 明 松 材 线虫 对 于 松 墨 天 牛 更 像 是 一 种 压力 因 
子 的 存在 ,而 不 是 一 种 外 来 侵 染 物 ,可 能 是 二 者 长 期 
















































































































































































协同 进化 的 过 程 中 建立 的 特殊 共生 关系 的 体现 。 

根据 KEGG 分 析 及 GO 分 析 可 看 出 , 当 松 材 线 
虫 进 入 天 牛气 管 后 ,表达 量 上 调 的 共 济 失调 蛋白 、 酌 
氨 酸 激酶 受 体 等 与 DNA 修复 和 细胞 损伤 修复 有 关 ， 
其 中 尤 以 PABPN1 蛋白 表达 量变 化 最 为 明显 。 而 
CREB 结合 蛋白 、 锌 指 和 蛋白 `.PABPN1 蛋白 .RNA Zi 
合 蛋 白 与 真 核 基因 表达 调控 以 及 RNA 合成 与 代谢 
有 关 等 方面 的 基因 表达 量 增加 ; 同时 作为 生殖 系 干 
细胞 调节 基因 的 核 小 核糖 核 蛋白 颗粒 SmE 则 与 生 
殖 活动 相关 (Ishigak et al., 2015; Dyson and Wright, 
2016; Banerjee et al., 2017) , 而 催化 活性 过 程 中 的 
上 调 的 FYVE 域 蛋白 则 与 激素 相应 和 应 对 胁迫 压力 
方面 有 潜在 作用 (Wang et al., 2000; Yao et al., 
2016), 。 同 时 不 难看 出 许多 上 调 的 基因 均 与 DNA 复 
制 \. 损 伤 修复 相关 ,如 Endonuclease V 5j DNA 的 损 
伤 修复 相关 ;Cystatin-C 基因 的 表达 产物 被 认为 对 半 
胱 氨 酸 和 蛋白酶 的 局 部 调节 有 着 十 分 重要 的 作用 ，; 
Putative helicase mov-10-B. 1 作为 一 种 可 能 的 RNA 
解 旋 酶 ,是 miRNA 介 导 的 基因 沉默 过 程 所 必需 的 ; 
而 Heat shock protein 75 kD, mitochondrial ( precursor) 
表达 形成 的 TNFR 及 其 下 游 激酶 则 直接 与 细胞 增殖 
调节 与 胁迫 响应 有 关 ( Copeland, 2012) 。 所 以 我 们 
推测 ,作为 松 材 线虫 的 传播 媒介 , 松 墨 天 牛 在 携带 松 
材 线虫 后 ,通过 增强 遗传 物质 与 细胞 组 织 修复 及 真 
核 基因 转录 调控 等 方面 的 基因 表达 ,来 提高 自身 对 
于 松 材 线虫 的 耐 受 性 ,抵消 松 材 线虫 寄生 所 带 来 的 
不 利 影响 ,重新 维持 了 自身 正常 生命 稳 态 ,从 而 为 松 
材 线虫 与 松 墨 天 牛 之 间 形 成 共生 关系 提供 了 可 能 。 

本 研究 从 转录 组 测序 方面 探究 了 携带 与 不 携带 
松 材 线虫 状态 下 的 松 墨 天 牛气 管 .表皮 及 脂肪 体 3 
个 组 织 间 表 达 基 因 的 变化 ,为 探究 二 者 之 间 的 共生 
关系 提供 了 基础 数据 支持 ,为 松 材 线虫 防治 工作 提 
供 了 新 的 参考 思路 。 后 续 我 们 将 对 重点 差异 表达 基 
的 功能 进行 分 析 , 更 全 面 深入 地 解析 松 墨 天 牛 在 
携带 松 材 线虫 后 所 发 生 的 生理 变化 的 分 子 机 制 ,为 
将 来 的 实际 应 用 打下 良好 基础 。 
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